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ABSTRACT 


This report considers empirical solutions to the first-passage 
or single-highest-peak (SHP) problem for a distributed elastic structure 
with rectangular geometry subjected to both stationary and a specific 
form of nonstationary random excitation. The structure is a flat, homo- 
geneous, uniform, square plate and the applied stationary excitation 
is white noise perfectly correlated in both space and time. The non- 
stationary excitation is a rectangular noise burst with a unity correla- 
tion in both space and time. The structure and the excitation are simu- 
lated electrically and peak response data are collected for (1) simply 
supported boundary conditions, at all edges and (2) rigidly clamped 
boundary conditions at all edges. These response data are used to 
establish probability curves yielding an estimate of the probability that 
the maximum response, for a finite time interval, remains below a pre- 
selected threshold level. 
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1. INTRODUCTION 


The SHP (single-highest-peak) problem generally refers to predict- 
ing the maximum response a mechanical system may achieve, within a 
finite time interval, when subjected to stationary, broad-band white noise 
excitation. It is noted that the maximum response need not necessarily 
correspond to a peak amplitude in the response time history. The solu- 
tion rests with determining ( | (3 j < (3^ , T) which is the probability that 
the absolute value of the response (3 remains below a threshold level (3^ 
within the time interval T . The term (3 is a dimensionless ratio of the 
maximum amplitude to RMS response of the system and T is the time 
interval over which the response is observed. The (3 response may be 
a positive maximum value p + , a negative maximum value p , or an 
absolute maximum value |p| and the RMS response is the system RMS 
response to stationary white noise excitation perfectly correlated in 
space and time. For this discussion, only solutions involving |p| 
values are shown. 

For a single degree-of-freedom system, the SHP problem has 
been examined in detail by various investigators using both simulation 
and analytical procedures as is discussed in References listed at the 
end of this report. Solutions acceptable for practical applications have 
been cleverly determined not only for stationary white noise excitation 
(References 1 and 7), but also for white noise applied as a step function 
(Reference 5) and for particular noise bursts as well (References 4 and 8). 
It is judicious to examine some of the established results for the vibra- 
tion of a single degree-of-freedom system before discussing the SHP 
problem for a distributed elastic structure. 
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2. RESULTS FOR A SINGLE DEGREE-OF-FREEDOM SYSTEM 


A linear single degree-of -freedom mechanical system (also called 
a linear mechanical oscillator) is shown as Figure 1 



Figure 1. A Single Degree-of-Freedom Mechanical System 


The dynamic behavior of this system is defined by the equation of motion 


mx + cx + kx = f(t) (1) 


where m corresponds to the mass, c the viscous damping coefficient, 
k the spring constant and f(t) a forcing function applied to the mass. 

In an alternate form, this equation may be written as 
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( 2 ) 
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and £, is the damping factor and oo the undamped natural frequency of 


n 


the system. For f(t) equal to a simple harmonic forcing function, 
say f(t) = f sin oot , the forced displacement response x(t) is 


^ H(co) r 

x(t) = - f sin oot 
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where the magnification factor H(oo) is defined in complex form by 


H(oo) = 
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For f(t) equal to an arbitrary deterministic function of time, the re- 
sponse may be calculated in the following ways 


x(t) 



H(w) F(w) e 1Wt dco 



h(T) f(t - t) dT 


( 6 ) 
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where F(co) is the Fourier Transform of f(t) and 1 i(t) is the system 
response to a unit impulse excitation. The solution obtained by inte- 
grating over co is a Fourier integral solution while the form given by 
integrating over t is a convolution integral solution. The Fourier 
integral gives the solution x(t) as a superposition of steady-state re- 
sponses and the convolution integral gives the solution as a super- 
position of free vibration responses. 

The output response of a lightly damped single degree- of-freedom 
system wherein f(t) is stationary broadband white noise excitation is 
shown as Figure 2. The undamped natural frequency of the mechanical 
system is assumed large compared with the half-power bandwidth AND 
the bandwidth of the input spectrum is assumed wide in comparison to 
the system half-power bandwidth and includes f . It is further implied 
the system has achieved stationarity in its response and we are ob- 
serving the response beginning at an arbitrary instant of time. 


Amplitude 



Time 


Figure 2. Response of a Lightly Damped Mechanical Oscillator to 
Broadband Stationary Random Excitation 
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This response appears as a sinusoid (at the frequency f ) with 

n 

a slowly varying random amplitude and random phase. The average 
number of zero crossings per second E [N^] is approximately 


E[ N o ] 
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where co and to are interpreted as the half-power point frequencies 

1 1 Z 

in radians per second. The mean square displacement response 1 is 


where 
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and is the magnitude of the input spectrum with units of lb /rad/sec. 

The term I (see page 80 of Reference 6 and Reference 11) is a dimension- 
n 

less coefficient ranging from zero to one with its exact value being de- 
pendent upon (1) the damping factor £ and (Z) the ratio of the cutoff 

frequency to of the input spectrum to co as shown in Figure 3. Simi- 
c n ^ 

larly, the mean square velocity response is 



ttG Q 
o 


Z m co 

n 


II 

n 


(9) 


5 


3. 


FREQUENCY RATIO 



Figure 3. Numerical Value of I 
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where the dimensionless coefficient II is shown as Figure 4. For 

high Q systems (£, < 0. 010) , both I and II are nearly zero for co / go 
° n n c n 

somewhat less than one; they rapidly approach unity as co^/co^ becomes 
slightly greater than one and converge to unity as co^/co^— ** oo . This be- 
havior implies the mechanical system acts as a highly selective band- 
pass filter for both displacement and velocity and admits frequency com- 
ponents principally near co . By selecting I and II as unity, upper 
^ 2 Z 

limits for d/ and ib « are calculated, 
x x 

The autocorrelation function for the displacement response is 
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where co ^ is the damped natural frequency of the system and related 

/ r 

to co^ as co^ = y 1 - £, . The magnitude of this function at t = 0 yields 

the mean square displacement response 
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which is identical to Eq. (8) with 1^ equal to unity. 


•V 

The mean square acceleration response may be developed in a similar 

manner and, in general, tends to be divergent as co /co >1 and becomes 

unbounded as co /co — * oo . 

c n 
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MAGNITUDE OF II 





The time lapse At for the amplitude of R^(t) to decay to — times its 
initial amplitude is given by 


At 
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and is called the correlation time or correlation interval (Reference 5). 
This time interval corresponds to the time lapse for the free vibration 
of a lightly damped single degree- of-freedom system to decay to —■ 
times an initial amplitude . 

For a white noise excitation, the probability density function of 
the response at any instant of time is Gaussian (Normal), The probability 
density of the envelope of the peaks is Rayleigh. The probability density 
of the individual peak amplitudes is, likewise, nearly Rayleigh and be- 
comes precisely Rayleigh when the damping becomes sufficiently small 
to define an infinitely narrow bandpass at f . 

The probability density of the response maxima during a finite ob- 
servation time is neither Gaussian nor Rayleigh and forms the crux of the 
SHP problem. A solution for this probability distribution function is 
given in Reference 10 as 


y\ 
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(|(3| < p T) — A 
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T > T 


cor 


( 12 ) 


where is a constant whose magnitude is dependent upon the initial con- 

ditions of the system, and t denotes the time lapse for the autocorre- 

cor 

lation function to decay to a negligible value. For most applications with 
(3 q > 2, A^ may be assumed equal to unity with trivial error. The 
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EXPONENTIAL CONSTANT 









term a is dependent upon the system damping Q and a preselected 

level of the p ratio. Plots of a versus Q are shown in Figure 5 
o o _j_ 

as families of curves in p^ for both p and |p| • This figure in con- 
junction with Eq. (12) provides a theoretical solution to the SHP prob- 
lem for a linear mechanical oscillator. 

Empirical solutions to this problem have been determined using 
techniques amenable to digital and analog computation. By means of an 

active electrical analog simulation (differential amplifier circuits), the 

* 

results shown as Figures 6, 7, 8 and 9 have been obtained. Figures 6 
and 7 are families of plots in Q of the dimensionless time parameter 
f T / Q versus | p| for P^( |p| < p^) = 0.50 and 0.95. The s e data may 
be compared directly with the results of Figure 5 by restating Eq. (12) in 
the form 


a. 
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ln[P M (|P| < P o >] 


2ir Q 
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f T 
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and evaluating this expression in the following manner 


These data are partially from Reference 4 and partially from inde- 
pendent analog studies carried out after the writing of Reference 4. Addi- 
tional data were obtained from a digital simulation study sponsored by 
Mr. E. H. Schell of AFFDL, Wright- Patters on Air Force Base, Ohio. 
These digital results tend to substantiate the shown analog data although 
the P results for the higher probability values were consistently lower 
(within 10% and generally not more than 5%). 


11 









a 

o 


In [0. 501 



a 

o 


In TO. 95] 



The data values for (3^ = 3 and p = 4 in Figure 6 and for (3 = 3 , 

(3 =4 and (3^ = 5 in Figure 7 agree closely with their equivalent values 

in Figure 5. It may thus be concluded that these analog results defini- 
tively support the theoretical solution of Eq. (12). 

Figure 8 is a set of plots in Q of | (3| versus f^T/Q for 
P M (|P| i P Q ) = 0.95 and -P^(|P| 5 P Q ) = 0.50. Only curves for Q = 20 
and Q = 50 are shown since these values are considered to be more 
typical of damping in physical structures, and thus, of most interest. 
These curves are exploded views of Figures 6 and 7 for f n T/Q < 5 and 
find application wherein the time interval T may be relatively short. 

The larger |p| values are associated with the lower damping ratios for 
all f^T/Q with the sharpest rate of increase noted for f^T / Q ranging 

from 0 to 3 . It is observed that 6 values of 3 for P_ _(| (3 I <6 ) = 0. 50 

^ r o M 1 1 — r o 

and p values of 4 for P^.(| p| p^) = 0. 95 are to be expected where 

f T/Q > 5. 
n — • 

Figure 9 shows |p| versus f^T/Q plots of the average maxima of 
|p| for Q = 20 and Q = 50 . These curves are very nearly the same as 

A, 

those for P^.(jp| < P Q ) = 0. 50 in Figure 8 with one exception. The data 
for Q = 20 appears to become numerically equivalent with that for Q = 50 
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Figure 9. Average Response Maxima of Single Degree-of-Freedom System 
to Stationary White Noise 




where f^T/Q > 4. From these data, the median and average values of 
lel for Q = 50 tend to become very nearly the same. It is suggested 
that a practical estimate of the average maximum response may be de- 
termined by using the|(3| versus f T/Q plot for Q = 50 . 

n 

The maximum response of a simple oscillator to shaped noise 
bursts also has been examined by an active analog simulation as dis- 
cussed in References 4 and 8. The noise burst excitation may be written 
as 


f(t) = E(t) n(t) 


(14) 


where E(t) is a specifically defined envelope function and n(t) is the 

output response of a random noise generator with a white noise spectrum 

of magnitude . Data from these studies are shown as Figures 10 

and 11 where is the dimensionless ratio of the maximum response 

to the RMS response and f t/Q is the dimensionless time ratio where r 

n 

is a time duration of the input burst of noise. The RMS response is the 
output response of the oscillator to broadband white noise wherein 
E(t) = 1 . By normalizing the (3 response in this way, time varying 
RMS considerations are avoided. These data suggest the noise burst 
response maxima become relatively the same as that for the stationary 
response when f^T/Q -1.5. Thus, plots such as these may be used to 
establish criteria for the conditions under which an oscillator achieves 
stationarity in its response. 
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3. RESPONSE MAXIMA OF DISTRIBUTED ELASTIC PLATES 


In this section, the SHP problem is considered wherein the struc- 
ture is a flat square plate and the applied force excitation is stationary 
white noise correlated perfectly in both space and time. The (3 re- 
sponse statistics are then examined via measurements taken at a point 
near the center of the plate for two types of boundary conditions, (1) 
simply supported at all edges and (Z) rigidly clamped at all edges. 

The equation of motion for a flat square plate is 


d V^( w ) + = f( x > y> t) 


(15) 


where f(x, y, t) is the applied force excitation per unit area, m the mass 
per unit area, w = w(x, y, t) = the lateral deflection from static equilibrium, 
v the Poisson's ratio and 


( > = 


d 2 ( ) 

dt 2 


D = 
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v 4 ( ) = 


ax 


8y 
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(16) 


The coefficient D is the flexural rigidity of the plate wherein E is 
Young's modulus and h the plate thickness. The notation ) 
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defines a spatial operator suitable to the rectangular geometry of the 
square plate where is the familiar dell operator for rectangular co- 
ordinates. 

The mean square displacement response of an arbitrary linear 
elastic structure to a uniform, homogeneous forcing field may be 
written as 



■j' 

H.(oo) H (go) Li (go) dco 
J k J k 


( 17 ) 


where the quantity L (go) 

jk 

spectral density function 


is expressed in terms 
S^(r, *r go) as 


of the spatial cross- 
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or in terms of the joint acceptance j (go) as 

jk 


S (7 cu) A 

= - A — z j ik (w) 
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The notation T refers to a spatial position vector, <(> (r) the jth normal 
mode of the system, coc the jth modal frequency, tT(co) the jth modal 
magnification factor defined as 
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and H (go) the complex conjugate of H (go) . The symbol A denotes the 

iC rC 

area of the plate and the quantity S^(r^, co) refers to a spectral density 

function of the applied force at the spatial location r* . If the position 

o 

r is selected to be a point on the structure where S.(r , oo) is a 

° 2 f ° 
maximum, j (go) will vary from zero to one. From Eqs. (18) and (19) , 
jk 

the joint acceptance is found to be 


j 2 k (co) = 
3 k 


STr , oj) A 
f o 


f A f\ 

J o J o J 


(rj <t> k (r') S f (r, r', to) dr dr - ' 


( 21 ) 


By assuming the j / k terms to be negligible in comparison to the 
j = k terms, the mean square displacement response becomes (for the 
conditions of this problem) 



or, 


in terms of the physically realizable spectrum , 
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Similarly, the mean square velocity response may be written as 
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In terms of the functions I and II shown as Figure 3 and 4, Eqs. (23) 

n n 

and (24) are of the form 
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The SHP problem for the plate structure is studied empirically 
by examining the (3 results from the electrical analog simulation de- 
picted schematically by Figure 12. The transformer whiffle-tree cir- 
cuit is used to create a spatial correlation of unity for the white noise 
excitation which is applied to the plate analog. As shown in this figure, 
the data of interest include both positive and negative response maxima 
(measured by the "peak" readers) and the RMS response. 
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Excitation Velocity 

Response 


Figure 12. Block Diagram of Analog Simulation Study 


The plate analog circuits are derived from energy concepts (Ref- 
erence 2) and may be categorized as mobility analogs. Such circuits 

consist of passive electrical components, are appropriate for nonuniform 

* 

physical and geometric properties, and appear topologically similar 
to the physical system. Mathematically, this type of analog corresponds 
to a finite -difference model whereas, mechanically, it may be interpreted 
as a form of lumped-parameter model. Background information about 
such structural analogs is found in Reference 9. 

The physical dimensions and physical properties of the uniform, 
homogeneous plate are defined as 


*T* 

For the more complicated structural configurations, a rather vivid 
imagination admittedly is required. 
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a = b = i = llin 


h = 0. 0625 in 
E = 10. 5 x 10 6 lb/in 2 
p = 0 . 1 0 lb / in^ 
v = 0. 3 


For simply supported edge conditions, the plate modal frequencies are 
given by 



(27) 


where rs denote the mode numbers (the number of half-sines| in the 

- 5 2 3 

x and y dimensions). With D = 235 lb-in and m = 1. 62 x 10 lb- sec /in , 
the fundamental modal frequency is f — 99 cps . For rigidly clamped 
edge conditions, the fundamental modal frequency is f^ —181 cps. 

The plate circuits and expressions for the component values used 
in this study are given in Reference 3. By scaling the analog circuits 
such that electrical time corresponds to real time (N = 1) and applying a uni- 
formly distributed sinusoidal load, the velocity to force frequency response 
4i 4i 

functions at x = — • and y = — ■ are obtained as shown in Figure 13. 

9 y 9 * 

These response functions indicate the magnitude of the circuit impedance 

at a particular point and are equivalent to mobility plots (magnitude only). 

The fundamental modal frequencies are noted to be approximately equal 

to the values previously calculated for both simply supported and rigidly 

clamped boundary conditions. Due to the symmetry of the loading, only 

even valued modes (i. e. where the mode shapes are symmetric about the 
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mid-spans of the plate) can be excited as is seen in these plots. The 
depth of the valleys between the modal frequencies is indicative of the 
modal coupling in the system. 

A measure of the damping in the first mode is given by 


_ 11 
Q 11 “ Af 


11 


(28) 


where Af^ is the half-power bandwidth centered at f^ . This expression 
yields — 33. 9 for the simply supported boundaries and ~ 39* 7 

for the rigidly clamped boundaries. Such damping is inherent in the com- 
puter (due to the fact that the circuit capacitors, inductors and trans- 
formers are not ideal, inter connecting wires have resistance, physical 
layout of the components produce parasitics, etc. ) and, although it can- 
not be precisely controlled, damping values in the first mode ranging from 
Q = 30 to Q = 50 are to be expected with proper scaling. 

In the simulation study, the plate is continuously excited by station- 
ary white noise correlated perfectly in both space and time and the velocity 
41 4 i 

response at x = and y = — — is observed for various time durations T 
such that the dimensionless ratio f^T/Q^ ranges from 0 to 5. 00. 

The magnitude of the white noise input is adjusted such that a response 
maxima of approximately six times the RMS response (6 (3^ values) could 
occur without distortion. For a specified T value, the response is 
sampled 100 times wherein both the positive and negative response max- 
ima are recorded for each sampling thus producing a total of ZOO data 
points. Typical response time histories (as observed over a finite sam- 
pling time interval T) are shown as Figures 14 and 15. It is to be noted 
the response does not appear to be the same as the stationary response of 
a linear oscillator to white noise, i. e. , as a type of harmonic motion with 
randomly varying amplitude and phase. Rather, higher frequency effects 
due to modal frequencies other than f^ are evident. 
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Stationary White Noise Excitation 



Velocity Response 


Figure 14. Response of a Square Plate to Stationary White Noise with a Spatial Correlation of Unity 
Simply Supported Boundaries 
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Stationary White Noise Excitation 


Figure 15. Response of a Square Plate to 
Rigidly Clamped Boundaries 




Velocity Response 


iarv White Nni.qp with a ^vnatial Correlation 






These individual maxima response data are next arranged in 
rank order, then divided by the RMS response to form histograms of 
Pj^.(P < P ) and P^.(P < P ) • A typical example of such data is 

shown as Figure 16. By adding the (3 and P data, a histogram of 

A 

yipl < P ) is formed as shown in Figure 17 and this is used (one 
such histogram for each sampling time T) to obtain individual prob- 
ability data points for |p|^ versus f^T/Q^i plots. It is understood 
that these histograms are estimates of the true maxima response prob- 
ability density functions for a specified value of f T/Q^ . Additional 
data should therefore be accumulated so that the variance of such data 
may be analyzed. 

Following such an empirical procedure just described, the |p| 
versus f^T/Q^ plots of Figures 18 and 19 are developed. Note that 
only the fundamental modal frequency f and its associated damping 
Q need be known for the structure. These curves display a be- 
havior similar to that for the single degr ee - of-f reedom system (com- 
pare with Figures 8 and 9) but have consistently higher values of |pl 
for a given probability. The curve for P^. <1 pi < P ) = 0. 50 is approxi- 
mately the same as that for the average [ p | and consequently is not in- 
cluded as a separate figure. 

For a given probability, the |pj_^ ratio increases indefinitely with 
time and theoretically |p| — ► oo as T oo , the latter being a truism 
of limited practical concern. The greatest rate of increase in |p| 

is seen to occur within 0 < f n T/Q^ < 3 . For f T/Q — 3, jpf 
values of approximately three are expected on the average and 95% of 
all amplitude maxima are expected to be equal to or less than P^ — 4 . 
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Figure 18. Maximum Response of a Square Plate to White Noise Correlated Perfectly 
in Space and Time 


i 



Figure 19. Maximum Response 
Perfectly in Space a 






Consider the maximum response characteristics of a plate sub- 
jected to a burst of white noise shaped in time by a well defined envelope 
function. It is required to develop a solution such that one may predict 
the expected single highest response maximum due to a noise burst of 
finite duration. This problem may be examined in much the same way 
as with the stationary excitation, that is by an analog simulation study 
wherein data for the response maxima are collected and arranged to 
form | (3 | plots. However, special attention must be given to (1) the 
RMS response used to form the (3 ratio and (2) the time duration com- 
mon to the dimensionless time parameter f T/Q . As with the 
previous problem, the (3 response is normalized using the RMS response 
of the system to broadband white noise correlated perfectly in both space 
and time wherein the magnitude of the stationary noise corresponds to 
the value where the envelope function is a maximum (E(t) = 1) . Al- 
though admittedly arbitrary, this procedure avoids consideration of a 
time varying RMS. The sampling time interval T is replaced by t 
which denotes the time duration of the noise burst. With these modifica- 
tions in definition, this nonstationary problem may now be plausibly 
attacked. 

With the system initially at rest, a rectangular noise burst is 
applied for t seconds typically producing the response as shown in Fig- 
ure 20. The system response is monitored from the onset of the ex- 
citation until the free vibration response has decayed to a negligible level 
and both the positive and negative maximum response values are re- 
corded. By sampling the system response one hundred times for each t 
setting and arranging the data in the same manner as with the stationary 
excitation, the |(3| plots of Figures 21 and 22 are constructed. Data 
points for simply supported boundaries are noted by • l s and data points 
for rigidly clamped boundaries are shown by A ! s . 
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Noise Burst Excitation 



Figure 20. Typical Response of a Square 
Correlated Perfectly in Space 








Figure 21. Maximum Response of a Square Plate to Rectangular Bursts of White Noise 
Correlated Perfectly in Space and Time 









Similar to the behavior for a single degree-of-freedom system, 
the response data are bounded by the curve depicting the stationary re- 
sponse of the plates. As the time duration of the noise burst becomes 
long relative to the fundamental modal period of the structure, the non- 
stationary excitation appears as a stationary white noise input and the 
noise burst response becomes essentially that of the stationary curve. 
Curves such as these may be used to define the conditions under which 
a multimode system excited by a noise burst achieves stationarity in its 
response. 

Of vast importance, however, is that the data for both stationary 
and nonstationary random excitation imply the I Pi sponse tends to be 
independent of the boundary conditions. Such a statement must be con- 
sidered as a qualified conjecture until more data are collected and ana- 
lyzed and/or theoretical work is advanced to support or reject this 
claim. In the interim, lei curves such as included in this report may 
be used discretely in engineering applications wherein a maxima re- 
sponse criteria for distributed structures is applicable. 
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4. CONCLUDING REMARKS 


1 


The SHP problem has been reviewed briefly for a single degree- 

of-freedom system wherein both analytical and empirical solutions have 

been discussed. More important, data have been presented forming an 

empirical solution to the SHP problem for a square plate with (1) simply 

supported boundaries and (2) rigidly clamped boundaries. The random 

excitation applied over the plate is (1) stationary white noise with a 

unity correlation in both space and time. and (2) nonstationary excitation 

shaped (in time) as rectangular noise bursts with a unity correlation in 

/s 

both space and time. A comparison of < P Q ) data for the single 

degree-of-freedom system and plate shows the IpI response maxima to 

be consistently higher for the plate. Of particular significance in the 

/\ 

plate data is the implication that the < P Q ) curves are inde- 

pendent of the boundary conditions. Such an implication has far-reaching 
practical significance and should be investigated thoroughly in the im- 
mediate future. In addition, the effect on the response maxima of mass 
loading and of bandlimited excitation with spatial and spectral character- 
istics more typical of service conditions would be indeed worthwhile 
investigations. 
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